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Abstract 



We report on a calculation of Higgs-boson contributions to the decays 
Bg — > and B — > Kl~^l~ {I = e, /i) which are governed by the effec- 

tive Hamiltonian describing h — > sl~^l~ . Compact formulae for the Wilson co- 
efficients are provided in the context of the type-ll two-Higgs-doublet model 
(2HDM) and supersymmetry (SUSY) with minimal flavour violation, focusing 
on the case of large tan (5. We derive, in a model-independent way, constraints 
on Higgs-boson-mediated interactions, using present experimental results on 
rare B decays including b sj, Bs — > Ai"*"/^", and B K^*^ ijl^ ^~ . In par- 
ticular, we assess the impact of possible scalar and pseudoscalar interactions 
transcending the standard model (SM) on the branching ratio of Bg Ai"*"/^" 
and the forward-backward (FB) asymmetry of fi~ in 5 — s- KjjL^ decay. 
The average FB asymmetry, which is unobservably small within the SM, and 
therefore a potentially valuable tool to search for new physics, is predicted 
to be no greater than 4% for a nominal branching ratio of about 6 x 10^^. 
Moreover, striking effects on the decay spectrum oi B ^ Kfi^n^ are already 
ruled out by experimental data on the Bg — > ix^ fi" branching fraction. In 
addition, we study the constraints on the parameter space of the 2IIDM and 
SUSY with minimal flavour violation. While the type-II 2HDM does not 
give any sizable contributions to the above decay modes, we flnd that SUSY 
contributions obeying the constraint on 6 — > 57 can signiflcantly affect the 
branching ratio of Bg — > ^"'"/x". We also comment on previous calculations 
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contained in the literature. 
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I. INTRODUCTION 



At the quark level, the decays Bg —>■ l^l^ and B Kl^l^, where I denotes either e 
or n, are generated by the short-distance effective Hamiltonian for b sl~^l~!^ Within the 
standard model (SM), the decay Bg l^l^ proceeds via penguin and box-type diagrams, 
and its branching ratio is expected to be highly suppressed. Likewise, the forward-backward 
(FB) asymmetry of the lepton in i? ^ Kl^l~ is exceedingly small. However, in models 
with an extended Higgs sector these observables may receive sizable contributions, and thus 
provide a good opportunity to look for new physics. The models to be considered are a 
type-II two-Higgs-doublet model (2HDM) and a supersymmetric extension of the SM with 
minimal flavour violation (see, e.g., Refs. - that is, we assume the Cabibbo-Kobaya- 

shi-Maskawa (CKM) matrix to be the only source of flavour mixing. An interesting feature 
of these models is that large values of tan/3, the ratio of the two vacuum expectation values 
of the neutral Higgs fields, may compensate for the inevitable suppression by the mass of 
the light leptons e or /i. 

The calculation of Higgs-boson exchange diagrams contributing to the h — > sl^l~ tran- 
sition has been the subject of many investigations PHTT|. As pointed out in Ref. 0, the 



results obtained in the context of the 2HDM disagree with each other. In view of this, 
we re-analyse the h sl^l~ transition, confining ourselves to the case of large tan/3 in 
the range 40 ^ tan/3 ^ 60. Our study extends previous analyses in several ways. We 
include, for example, other rare B decays in addition to 6 57 to constrain possible scalar 
and pseudoscalar interactions outside the SM. We also assess their contributions to various 
observables in 6 ^ sl^l~ transitions such as i? — > K^*^l~^l~. 

So far all experimental results yield only upper bounds on the decay modes governed by 
b —>■ sl^l^ . The best upper limits at present come from processes with muons in the final 
state, and we therefore concentrate on the fi~^fi~ mode. Specifically, we address the viability 
of the short- distance coefficients in the presence of scalar and pseudoscalar interactions with 
the measured 6 — 57 rate and the experimental bound B{Bs /^^/^~) < 2.6 x 10~^ 
(95% C.L.) [|12|, as well as with the restrictions imposed by the upper limits on B ^ 

The outline of the paper is as follows. The effective Hamiltonian describing the quark 
transition b sl^l~ in the presence of non-standard Higgs bosons is reviewed in Sec. 0. 
In Sec. |Tl|, we discuss the hadronic matrix elements required for the decays Bg —>■ and 
B — »• Kl^l~ . The corresponding angular distributions and decay spectra are presented in 
Sec. |V[ Section |V| is devoted to the calculation of the Higgs-boson diagrams in a general 
i?^ gauge, and also contains a brief description of our renormalization procedure. Readers 
who are not particularly interested in the details of the computation, can skip this part and 
proceed to the discussion of our results, which are obtained in the framework of the type-II 
2HDM and supersymmetry (SUSY) with minimal flavour violation. Attention is focused on 
the interesting case of large tan /3, and a comparison is made with the results of previous 
studies. In Sec. |V|, we derive model-independent upper bounds on scalar and pseudoscalar 
interactions, and explore their implications for the branching fraction of -B ^ K^'^fi~, 



^We adopt a convention where B = (db) or [ub) and Bg = (sb). 
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as well as the corresponding FB asymmetry, using presently available data on the decays 
fi'^fj, , and Bg — >■ . As an application, we investigate the constraints 

on the parameter space in the aforementioned extensions of the SM. Finally, in Sec. |V11| , we 
present some concluding remarks. 



II. EFFECTIVE HAMILTONIAN FOR b sl+l' 



The starting point of our analysis is the effective Hamiltonian describing b sl^l : 

AG r ^° 

* ^ i=l 

+ c's{ii)0's{ii) + dp{i,)0'M^, (2.1) 

where cf\fi) and 0^'\fi) are the Wilson coefficients and local operators respectively. In 
writing Eq. (|2.1| ), we have used the unitarity of the CKM matrix and omitted terms pro- 
portional to V^^V*jVf.^V* ~ 0(10-2). Taking the limit cfp 0, we recover the effective 
Hamiltonian of the SM |T|,|l|]. 

The evolution of the short-distance coefficients evaluated at the matching scale nw = 
down to the low-energy scale at fXb = ml°^'^ can be performed using renormalization group 
equations (see, e.g., Ref. [^). The operator basis is given by 



Ol = {Sal^.PLCp){cp^^PLha), 
q=u,d,s,c,b 

0, = {s^^,PLbf,) Yl ((IP^'^PLQa), 

q=u,d,s,c,b 

0, = {s^^^PLba) Yl i'iP^'^Pm)^ 

q=u,d,s,c,b 

0, = {-s^l^PLbp) Y (Ip^'PrI-), 

q=u,d,s,c,b 

^7 = T^rnb{s^a^uPRba)F^''', 

^8 = ^MSaT^f,^^.PRbf,)G'^^\ 
2 



IGvr 



0, = ^is^rPLb^)ih^l), 

0,0 = ^,{Sal^PLbo){h,l^l), 
Os = Y^^b{SaPRba){ll), 

_ 

Op = -r-^rnb{saPRba)ih5l), 



O' 



167r2 



(2.2) 



a, /9 being colour indices, a labels the SU(3) generators, and Pl,r = (1 =F75)/2. Notice 
that we have dropped the rUg corrections to O-i and Os while retaining the primed operators 
O'g and O'p. Further discussion on this point will be given in Sec. 0. In general, there are 
additional operators such as {sa ^yPL^iih){la^^^ PlrI) which, as we will argue later, do not 
contribute to the decay Bs l^l^ but show up in the process B Kl^l~ . However, these 
operators are not expected to contribute significantly |TB], and we shall neglect them in our 
subsequent discussion. 



III. HADRONIC MATRIX ELEMENTS 



A. B ^ Kl+l- 



The hadronic matrix elements responsible for the exclusive decay B —>■ Kl^l are con- 



veniently defined as p!7|-p!9 



/\/f2 _ /If 2 

mk)\-s^,b\B{p)) = {2p - q),U{q') + \, ^ g.[/o(g^) - U{q% 



(3.1) 



{K{k)\-sia,,q^h\B{p)) = -[(2p - q)^q' - (M| - M|)g^] 



fAq') 



Mb + Mk' 



(3.2) 



where = {p — k)'^ is the four- momentum transferred to the dilepton system.Q Further, 
employing the equation of motion for s and b quarks, we obtain, from Eq. (|3.1| ), 



{K{k)\-sb\B{p)) = ^ ^fo{q') 

THb - rus 



(3.3) 



In the following we shall adopt the form factors of Ref . |I8| , which uses light cone sum rule 
results. 



B. Bs l+l' 



The relevant matrix elements are characterized by the decay constant of the pseudoscalar 
meson Bg, which is defined by the axial vector current matrix element 



(0 1 57^75^1 5s(p)) = ^P^^fBs 



(3.4) 



^Note that {K\sj^{l - -f5)b\B) = {K\sj^{l + j5)b\B) = {K\s-f^b\B). 
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while the matrix element of the vector current vanishes. (For the decay constant we will 



take the value 210 ± 30 MeV |20|.) Contracting both sides in Eq. ( p.4|) with p^^ and using 
the equation of motion gives 

{0\sj,b\B,{p)) = -ifs^ . (3.5) 

nib + nis 

An important point to note is that the matrix element in Eq. ( p.4|) vanishes when contracted 
with the leptonic vector current /7^/ as it is proportional to = p'^^ which is the only 

vector that can be constructed. In addition, the matrix element {Q\sa ^yh\B s{p)) must vanish 
since it is not possible to construct a combination made up of that is antisymmetric with 
respect to the index interchange fi ^ u. Consequently, the operators O7 and Og do not 
contribute to the decay Bs —>■ l^l^, which is then governed by Oio and O^^p defined in 
Eq. (pD. 



IV. DIFFERENTIAL DECAY DISTRIBUTIONS 



Using Eq. ( p.l| ) together with Eqs. ( |3.1j )-( |37^ ), the matrix element for the just-mentioned 
decay modes can be written in the form 



Fsll + Fph^l + + FAP^h^^l^l 



(4.1) 



p^^ being the four- momentum of the initial B meson, and the Fj's are functions of Lorentz- 
invariant quantities. It should be emphasized that the form factors Fs and Fp must vanish 
when = because of chiral symmetry, hence Fs^p cx m^. Nevertheless, as will be elab- 
orated below, large values of tan (3 may compensate for the suppression by the electron or 
muon mass in certain extensions of the SM. 

Squaring the matrix element and summing over lepton spins, we find the result 

\M\'' oc (s - 4m2)|F5p + s|Fp|2 + [4(p ■ pi-){p ■ p,+ ) - M^^ldF^P + |Fyp) + 4m2MV^p 
+ Ami[p ■ {pi+ - pi-)Re {FsF*y) + {p ■ g)Re (FpF^)], (4.2) 

where s = q^, q = pi+ + pi~, and M refers to the mass of the decaying B meson. 

A. B ^ Kl+l- 

Let us start with the decay B — Kl^l^, where we will employ the definitions 

X{a,b,c) =a^ + b'^ + -2{ab + bc + ac), pi = ^l- Amf/s. (4.4) 

Furthermore, we define 6 as the angle between the three-momentum vectors p^- and Ps in 
the dilepton centre-of-mass system. The two-dimensional spectrum is then given by 
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X [1 - /3f cos^ + |Fy n + Am^iMl\FA\^ + 2m,[A^/2(M|, M^, s)ARe cos0 

+ (M2 -Mi + s)Re(FpFi)]j, (4.5) 



(4.6) 



with s and cos 6' bounded by 

Amf ^ s ^ (Ms - MkY, -1 ^ cos^ ^ 1. 
A quantity of particular interest is the forward-backward asymmetry 

dV 

I dcosd- — ; — / dcosU— 

Ir. d.M rxcoR H I 1 d.M 









.1 dsdcosO 



dcos6- 



ds dcos 9 



•° dT 
+ I dcosO- 



(4.7) 



ds dcos 9 



which is given by 



^fb(s) 



dV /ds 



(4.8) 



where the dilepton invariant mass spectrum, dT/ds, can be obtained by integrating the 
distribution in Eq. (|4.5|) with respect to cos^. Exphcitly, we find 



^^^lE^lI^Hn = 2AV2(M|, Ml s)(3L{(3f\Fs\' + \Fp\^) + iA(M|, Mf, .) 
2ml 



X 1 + 



Fa\' + \Fvn + AmfMj,\FAY + 2m,(M^ - Mj, + s)Re {FpFl)] \. (4.9) 



The Lorentz-invariant functions in the above formulae depend on the Wilson coeffi- 
cients as well as the s-dependent form factors introduced in the preceding sections, namely, 



Fs = \{Ml-Ml)h{s) 



rrib - rus 



(4.10) 



— A/f2 

Fp = -miCro{ U{s) - ^fo{s) - U{s)] 



CpTTli, + dpTTls 



rrib - rris 



(4.11) 



F^ = cio/+(s), F, 



V 



cfU{s) + 2cfmb 



Ms) 



Mb + Mk 



(4.12) 



It should be noted that within the SM the Wilson coefficients of scalar and pseudoscalar 
operators are invariably suppressed by mimb^s/M^, leading to c^p ~ 0, and so the FB 
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asymmetry vanishes. This observable is therefore particularly useful for testing non-SM 
physics and its measurement could provide vital information on an extended Higgs sector. 

The analytic expressions for the remaining coefficients appearing in Eqs. ( ^4.10| )-( p:.12| ) 
may be found in Refs. |1^,|15|. Within the SM, they are estimated to beQ 



-0.310, 



„efr 



C9 + y(s), C9= 4.138, Clo 



-4.221, 



(4.13) 



where the function Y{s) denotes the contributions from the one-loop matrix elements of the 
four-quark operators Oi-Oq (see Appendix ^ . 

Finally, we give here the SM prediction of the non-resonant branching fraction for the 
decay into a fj^'^fi" pair, the result being 



B{B ^ Kfi+fi^) = (5.8 ± 1.2) X 10" 



(4.14) 



where the error is due to the uncertainty in the hadronic form factors, which is the major 
source of uncertainty in the branching ratio. We do not address here the issue of resonances 
such as J/ip, which originate from real cc intermediate states. For theoretical discussions 
of these contributions and the various approaches proposed in the literature, the reader is 



referred to Refs. [18, 22 . 



B. B, 



l+l- 



Our results for the matrix element squared [Eq. ( [4. 21 )] are immediately adaptable to the 
process Bg — > Using p = pi+ + pi-, we obtain the branching ratio 



B(5,^/+r) 



167r3 




1 - 



4mf 



\Fs\' + \Fp + 2miFA\ 



(4.15) 



The factor mi in front of Fa reflects the fact that within the SM the decays Bg e^e~ or 
/i"*"/!" are helicity suppressed due to angular momentum conservation; indeed, since the B 
meson is spinless, both and /~ must have the same helicity. 

The scalar, pseudoscalar, and axial vector form factors are given by {i = S, P) 



F, = --MIJb^ 



drrib 



nib + ms 



--/s.cio. 



(4.16) 



Throughout the present paper we use the leading-order result for the Wilson coefficient 
cio in order to be consistent with the precision of the calculation that will be presented in 
Sec. This is different from Refs. P,|ll|| where the next-to-leading-order result for the SM 
contribution has been used. 



^We use a running top-quark mass of rrit = mt{mt) = 166 it 5 GeV, corresponding to m 



pole 



174.3 ±5.1 GeV ^\ 
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For completeness, let us record the SM branching ratio for the dimuon final state: 

B {Bs = (3.1 ± 1.4) X 10"^ (4.17) 

where we have used the value | VJsl = 0.04 ± 0.002 along with the aforementioned ranges for 
fB^,mt. We emphasize that the error given is dominated by the uncertainty on the B meson 
decay constant fs^ ■ Before moving on to the computation of Higgs-boson exchange diagrams 
that contribute to the form factors F^, we briefiy recall the experimental constraints relevant 
to our analysis. 



C. Experimental constraints 

To date, the most stringent bounds on the magnitude of the previously discussed short- 
distance coefficients come from the Collider Detector at Fermilab (CDF) ||13|| : 

B (5° K*°fi+fi-) < 4.0 X 10"^ (90% C.L.), (4.18) 

which should be compared with the branching fraction of about 2 x 10^^ predicted by the 
SM. Also, from the absence of any signal from the process B^ — > K^fi^fj,', the 90% C.L. 
limit 

B{B+ ^ 12+ fi-) < 5.2x10-^ (4.19) 



has been derived |T3|, which is an order of magnitude away from the SM prediction of about 
6 X 10"^. 



The measurement of the inclusive branching ratio B {B ^ -^s7) yields the result ||23| 

2.0 X 10'^ <B{B ^ X,7) < 4.5 x 10^^ (95% C.L.), (4.20) 

which places limits on the absolute value of cf^. In what follows it is more convenient to 
define the ratio R-/ = Cj^/c^^' Using the leading-order expression for B{B ^ -^s7) from 
Ref. [53, we calculate the bound to be 



0.88 < |i?7| < 1.32. (4.21) 
A search for the decay Bg — > fi+^~ has been made by CDF, leading to the result |T^ 



B {B, /iV") < 2.6 X 10~*^ (95% C.L.). (4.22) 

This in turn translates, via Eq. (|4.15|) , into an upper limit on the strength of scalar and 
pseudoscalar interactions, as we shall discuss. 

We conclude this section with a few remarks on the B mode. Experimental search leads 
to a 95% C.L. (90% C.L.) upper limit oi B {B ^ < 8.6 x 10"^ ^ (6.1 x 10"^ [0), 

which is several orders of magnitude above the SM expectation of O(10~^°) We stress 
that if fiavour violation is due solely to the CKM matrix, the subject of the present paper, 
the B decay is suppressed relative to the Bs decay by a factor |Vtd/VtsP ~ O(10~^); however, 
this suppression does not pertain to models with a new fiavour structure. 
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V. HIGGS-BOSON CONTRIBUTIONS TO 6 ^ sl+l' 



We now turn our attention to the computation of Higgs-boson contributions to the Wilson 
coefficients of the scalar and pseudoscalar operators in the b —>■ sl^l~ transition, within the 
context of the 2HDM and the minimal supersymmetric standard model.0 

As anticipated at the outset of this paper, we evaluate the diagrams in the gauge, 
which provides a check on the gauge invariance of our calculation. We use the Feynman 
rules of Ref. pOl] and focus on the large tan/3 scenario, that is, 40 ^ tan/3 ^ 60. 



A. Two-Higgs-doublet model 

We compute the Higgs-boson exchange diagrams in the framework of a 2HDM where the 
up-type quarks couple to one Higgs doublet while the down-type quarks couple to the other 
Higgs doublet (usually referred to as model II), which occurs, for instance, in supersymmetry. 
We will use the SUSY constraints on the parameters A, appearing in the Higgs potential 
(see, e.g., Ref. 0). We defer the discussion of the more general 2HDM with Ai = A2, as well 
as the comparison with results presented in the literature, to the end of the section. 

The relevant Feynman diagrams for h sl^l~ are depicted in Fig. |l|, where and 
h^,H^ are the CP-odd and CP-even Higgs bosons respectively, represents the charged 
Higgs bosons, and are the would-be-Goldstone bosons. Before stating the results 

for the scalar and pseudoscalar Wilson coefficients, we pause to outline our renormalization 
procedure. 



1. Remarks on the renormalization procedure and the renormalization group evolution 

Writing the interactions of Higgs bosons with down-type quarks appearing in the 'bare' 
Lagrangian of the 2HDM and the minimal supersymmetric standard model (MSSM) in terms 
of renormalized quantities, we obtain the counterterms necessary to renormalize the theory, 
namely, in the one-loop approximation. 




ie {cos a; — sin a} 
2M]Y sin 9\Y cos /3 



2 



+ 



R 



(5.1) 



The Higgs-boson contributions to , Cg , cio in the massless lepton approximation can be found 
in Refs. |27-29|. For a non-zero lepton mass, there are also box diagrams with charged Higgs 
bosons which, at large tan/3, contribute only to the helicity-flipped operators Og ~ (s7^Pr6)(/7^/) 



and O'lQ ~ {s"l^PRb){l^f_i'y^l) [cf. Eq. (2.2)]; however, their contribution is negUgible for / = e or /i. 
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e {sin/3; — cos/3} 
2Mw sin 9w cos /3 



{^0; G°} - 



R 



rris 
2 



(5.2) 



a being the mixing angle in the CP-even Higgs sector, and 9w is the Weinberg angle. The 
quark field renormalization constants are defined through 



hL 



2 



bare 









sl 




sr 


hi 




{ bR 



bare 




(5.3) 



where U is the unit matrix. They can be determined from the two-point function of the sb 
vertex, which is given by 



2 



Pr 



6Zi + 6Z 



nis^Zh + mi,5Z, 



sb 



Pr 



(5.4) 



We have chosen an on-shell renormalization prescription in which the finite parts of the 
field renormalization constants are fixed by the requirement that the fiavour-changing b — *■ 
s vertex vanish for external on-shell fields.^ We have checked that our approach, based 
solely upon one-particle- irreducible diagrams, yields a gauge- independent result, which is in 
agreement with Ref. 0. 

Finally, a few remarks are in order regarding the renormalization group evolution. Re- 
ferring to Eq. ( |2.2| ), the masses of the light quarks, and m^, appear in the scalar and 
pseudoscalar operators rather than in the corresponding Wilson coefficients, and hence must 
be evaluated at the low-energy scale. The anomalous dimensions of the light quark masses 
and the scalar as well as pseudoscalar quark currents cancel, and so the anomalous dimension 
of the above-mentioned operators vanishes (see also the discussion in Ref. [PT|). Another 
method commonly used in the literature (see, e.g., Refs. ||10|,|ll|) is to absorb the light quark 



masses into the Wilson coefficients, which are determined at the high scale. In this case, 
the scalar and pseudoscalar operators have a non-vanishing anomalous dimension, and the 
values of the corresponding Wilson coefficients at the low scale are obtained by means of the 
renormalization group evolution. It is important to stress that both methods are equivalent. 



Since our analysis is performed at leading order in the electroweak couplings, the only quantity 
that has to be renormalized is the wave function. 



12 



2. Results for scalar and pseudoscalar couplings 



Retaining only leading terms in tan /?, our results can be summarized as follows: 

w _ „box _ m;tan^/9 

cs --Cp --J^-^BixH^^x,), (5.5) 

^peng _ ^peng _ 

^count ^ _^co.nt ^ _ , X.) , (5.7) 

4M^ sm fc'w' 

where the superscripts denote the box-diagram, penguin, and counterterm contributions 
respectively, and Xi = mf/M^. The functions B and C are listed in Appendix In 
deriving Eqs. (p.5|)-(pT7D, we have used the relation 

where Mfjo , M^o are the tree- level masses of the CP-even Higgs bosons. Note that we 
have chosen tan/3 and the charged Higgs-boson mass mfj± as the two free parameters in 
this SUSY-inspired scenario. Turning to the coefficients of the helicity-flipped operators 
c'g p, they are also proportional to tan^/3 but their contribution to the decay amplitude is 
suppressed by a factor of ms/m^ compared to cs^p, and hence can be neglected. 
Summing all contributions results in 



mi tan^ (3 In r ''^'h± 
M^sin^ ^ly 1 -r' ^ ~ ~mf 

(We will compare our findings with other recent calculations below.) 



B. SUSY with minimal flavour violation 

Since we consider a scenario with minimal fiavour violation, i.e. we assume fiavour- 
diagonal sfermion mass matrices, the contributing SUSY diagrams, in addition to those in 
Fig. in consist only of the two chargino states (see Fig. |^).0 It is convenient for the subsequent 
discussion to define the mass ratios 

2 2 TTh^ 

•^ki 9 ; Vai 9 ; ^ij 9 ) lO.iUI 

ml± ml± ml± 

Xi Xi Xj 

with i>k, Ua, and xf denoting sneutrinos, up- type squarks, and charginos. In terms of these 
variables and recalling Eq. (|5.8|), we obtain 



^The gluino and neutralino contributions are deferred to another publication |31]. 
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FIG. 2. SUSY contributions to the process b sl~^l , as described in Sec. with indices 
a = 1,2, ... ,6; i,j = 1,2; k = 1,2,3. The counterterms for the penguin diagrams are given in 
Eqs. (O) and (O). 



^ i,j=l a=l k,m,n=l xf 
X (^yaiU*2Vil±^Ui2Vji^Di{Xki,yai,Zji)'^, (5.11) 



„peng 
^S,P 



, 2^ij)5(2b^k7n 



{Mu)kk 
\/2m-± 



- ^TTO?^ ^ ^ ^ «3(...)«.(r-)„.„rL,J. (5.13, 



where 



'WV"-H± ^'-"Wl i=l a=l m,n=l 



Mu = diag(m„, m^, m^), (5.14) 



= ^ ^ . [V2MwVa{r''^\a - {Mu)nnV.2{r''''%a]Xmn, (5.15) 

2v2 sm"' 9w 

with the ratio of CKM factors Xmn = ymb^ns/^tb^Zy ^^e functions -Di,2,3 are hsted in 
Appendix 0. In writing the above formulae, we have used the unitarity of the CKM matrix 
and the squark mixing matrices (for definitions see Appendix 0). Note that the chargino 
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contributions vanish when all the scalar masses are degenerate, reflecting the unitarity of 
the mixing matrices. Another noticeable feature is that the leading term in tan/3 comes 
from the counterterm diagrams. 

Turning to the Wilson coefficients of the helicity-fiipped scalar and pseudoscalar opera- 
tors entering the operator basis in Eq. (|2.2|), we obtain 



/count _ r count /\ ^ \* \]* /count _ r count /\ \* \]* /p- i (-\ 

Cfi" — [Cs [^mn ^nm)\ 5 Cp — -[Cp {Aran ^nm)\ 5 1,0. iOJ 

which are 0(tan^/5), while the remaining coefficients are proportional to tan^/3. Recall that 
the contribution of the Wilson coefficients in Eq. (|5.16|) to the decay amplitude is propor- 
tional to msia,v? (3 and hence comparable in size to the mf,tan^/5 contributions, Eqs. ( p.9p , 
( p. 11 ), ( 5.12|) . On the other hand, this contribution is negligible when compared with the 



leading term m;, tan (3 and so will be omitted in what follows. 



C. Remarks on previous results 

We conclude this section by comparing our results with previous calculations in the 
literature |]9|-|TI|. To this end, we consider the MSSM as well as a general 2HDM with 
Ai = A2 for the coupling constants appearing in the Higgs potential . We discuss the two 
scenarios in turn. 

(a) Working in the framework of the MSSM, the Higgs sector is equivalent to the one of 
the 2HDM with SUSY constraints 0. The results of Huang et al. can be checked by 
exploiting the tree-level relation in Eq. ( |5.8| ). Reducing their expressions for the one- loop 
functions to the compact formulae given above, and after correcting numerous typographical 
errors, we agree with their results. Note that our anomalous dimension is equal to zero, due 
to the running 6-quark mass entering the definition of the operators [see Eq. (|2.2| )]. In order 
to compare our findings with those obtained by Chankowski and Slawianowska [jlT|, we 



specialize to the case M2 ^ so that m-± ~ |/i| and m-± \M2\. In this approximation, 
and retaining only contributions of the lighter chargino and the scalar top quark, our results 



are in agreement with Eqs. (33) and (34) of Ref. |Tl]. 

(b) In the context of the general CP-conserving 2HDM with the constraint Ai = A2, the 
set of free parameters consists of M/jO, M^q, M^o, mH±, as well as the mixing angles a and 
The necessary Feynman rules are listed in Ref. apart from the model-dependent trilinear 
Higgs couplings gn+H-hO and Qh+h-h^i which can be found, for example, in Ref. [^]. (We 
have rederived these couplings confirming the result given there.) In the limit of large tan (3, 
these couplings reduce to 

QH+H-ho — T-r- sina cos^ a^M'^o — M^o) tan/3[l + 0(cot /?)], (5-17) 

9H+H-H<i ^ ^sin^acosa(M^o - M^o) tan/3[l + 0(cot /?)]. (5.18) 

Our results for the trilinear couplings disagree with those presented in Eq. (27) of Ref. [Q. 
In addition, we would like to stress that within the general 2HDM the mixing angles a and 
f3 are independent parameters, contrary to the statement made in that work. (This has 
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also been pointed out in Ref. IT^].) Taking into account the Feynman diagrams due to the 
trihnear Higgs couphngs in Eqs. ( ^.17|) and (|5.18|) , we obtain 



cs 



mi tan^ (3 
AM^ sin^ Ow 



Inr ,(M' 
+ sin2(2a'^ ' 



HO 



Ml.] 



1 — r 



1 — r + In r 



mt 



(5.19) 



and Cp as given in Eq. (|5.9| ). The first term is in agreement with the calculation of Logan 
and Nierste and with the result obtained by Huang et al. ||]_ 
term, it is absent in the expression given in Ref. [H and differs from that of Ref. 



As for the a-dependent 
|1|. We 

caution that in models such as the general 2HDM with a complicated parameter space the 
subleading terms in tan (3 might be of the same order of magnitude as the leading ones for 
certain values of the parameters. 



VI. IMPLICATIONS FOR THE DECAYS Bs n+ fJ-' AND B Kfi+fj,- 

In this section we explore the consequences of the current upper limits on rare B decays 
discussed in Sec. ^for scalar and pseudoscalar interactions. In the quantitative analysis, we 
use rUb = rnb{mi,) = 4.4 GeV and neglect terms of order mg/mf,, which is certainly sufficient 
for our purposes. 



A. Model-independent analysis 

We start by analysing the constraints on scalar and pseudoscalar interactions, as well 
as on i?-physics observables, in a model-independent manner. To this end, we define the 
following dimensionless quantities 

Ri = ojcf, Rs,P = mbCs,P, (6.1) 

with cf^ as in Eq. ( [4.13| ). For the purposes of this analysis, we assume that there are no 
additional CP phases, besides the single CKM phase, so that the -R's in Eq. ( |6.1|) are real 
(remembering that we omit terms proportional to V^jjV*g/V^fjV*g -C 1). 



1. Bounds on scalar and pseudoscalar couplings 

The most severe constraints on scalar and pseudoscalar interactions arise, as we will 
argue shortly, from the upper bound on the Bg A*^/^" branching fraction, Eq. ( [4. 221) , 
which maps out an allowed region in the {Rs,Rp) plane. This is illustrated in Fig. ||(a), 
where we have chosen Jb^ = 210 ± 30 MeV and assumed i?io = 1 (i.e., the SM value for 
Cio)- We note that the allowed region in the {Rs, Rp) plane is fairly insensitive to the range 
—2 < RiQ < 2 implied by the present experimental bound on B {B ^ K* ji^ jj,^). This can 
be easily understood from Eq. ( |4.15|) , where the contribution of -Rio, or equivalently of cio, to 



the branching ratio is helicity suppressed. It is important to emphasize that the maximum 
allowed contribution of scalar and pseudoscalar operators to the B K*fj,~^fj,~ branching 
fraction is consistent with the experimental upper limit, Eq. ( [4.181 ). As will become clear, 
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FIG. 3. Constraints on the coefficients Rs and Rp as determined from the upper hmit on 
B (Bs fi^fi~). (a) For Riq = 1 and fs^ = 210 it 30 MeV. The dark region indicates the allowed 
region, (b) For the central value fs^ = 210 MeV and two choices of Riq: —2 (solid line) and 2 
(dashed line), as described in the text. 

the new-physics contribution due to Higgs-mediated interactions does not significantly alter 
the maximum allowed values of Rg and -Rio-Q 

Taking Rw = ±2 and fs, = 210 MeV, we infer from Fig. |(b) the interval -4 < Rs,p < 4 
for scalar and pseudoscalar coupUngs. 

2. Branching ratio and FB asymmetry in B ^ Kfi'^fj," 

We now assess the implications of scalar and pseudoscalar interactions for the branching 
ratio and forward-backward asymmetry. For the paper to be self-contained, we also provide 
the analytic expression for the B — > K*n~^^~ branching fraction. If we keep in mind that 
the Wilson coefficients are real, we obtain 

B{B ^ Kfi+fx-) = [0.042 + 2.846i??o + 2.730i?^ + 0.043i?^ + 0.181i?| + 0.182i?^ 

- 0.l32RpRio - O.GSGRjRq + 0.522i?9 - OMbRj] x 10"^ (6.2) 

B{B ^ = [0-015 + 0.922i?^o + O.SMg + 0.212i?? + 0.014/?| + 0.014i?| 

- 0.015i?pi?io - O.mRjRg + 0.177i?9 - 0.046i?7] x 10"^ (6.3) 



^The bounds on Rg also depend on the sign and magnitude of Ry. For details, see Figs. 9 and 10 
in Ref. and Fig. | below. 
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MM 



for the choice 



TABLE I. Branching ratios of the decays Bs — > n^^f^ and B K^*^ 
{R7, Rg, Rio) = (—1.2, 1.1,0.8) together with the present experimental upper hmits (see Sec. IV). 



{Rs,Rp) 



B{Bs 



B{B ^ KpL+n-) 



B{B ^ K*^i+ii- 



(0,0) 


1.8 X 10-9 


6.8 X 10- 


-7 


2.8 X 10-6 


(2,2) 


1.1 X 10"^ 


8.0 X 10" 


-7 


2.8 X 10-6 


(3,-2) 


2.1 X 10-6 


9.3 X 10- 


-7 


2.8 X 10-6 


(-4,0) 


2.5 X 10^6 


9.7 X 10- 


-7 


2.8 X 10-6 


Exptl hmits 


< 2.6 X 10-6 (95% C.L.) 


< 5.2 X 10-6 (9 


0% C.L.) 


< 4.0 X 10-6 (90% C.L.) 



where the -R's are defined in Eq. ( |6.1j ). We note that the limits on Rs and Rp from the 
upper bound on B {B Kjji'^jj,^) [Eq. ( [4.19| )] are numerically less stringent than those 
derived previously from Bg jj,^ ji^ . Consequently, |-Rs,p| ~ 4 is essentially the largest 
value that is possible. Some representative results for the branching ratios of Bg fi~^fi~ 
and B K^*^ fi^fi are summarized in Table |I[ As can be seen, the branching ratio of 
B — K*fi^fi^ decay is essentially unaffected by the presence of Higgs-mediated interactions 
since the contributions of Rs,p in Eq. ( |6.3|) are suppressed compared to those of -R9,io. 
Furthermore, it is important to emphasize that large effects of scalar and pseudoscalar 
interactions on the B —>■ Kfi~^fi~ decay rate are already excluded by the upper limit on 
the Bg — >■ branching fraction. (This constraint has not been taken into account in 

the analysis of Ref. [0.) As for the asymmetry, our main interest is in the average FB 
asymmetry, which can be obtained from the expression in Eq. ( [4.8|) by integrating numerator 
and denominator separately over the dilepton invariant mass, leading to 

(Afb) = i?s(0.512 + 5.424i?9 - 0.681i?7) ( j ■ (6.4) 

To gain a maximum FB asymmetry, we fix Rs = —4 and Rp = allowed by current 
experimental data on B{Bs /i+/i-), together with the SM value of Rj = 1. Referring to 
Fig. ^, it is evident that the average FB asymmetry in 5 — > Kfi'^fi" decay amounts to ±4% 
at most, the actual value depending on Rg and Riq. We emphasize that some of the values 
of {Rq, Rio), while respecting the upper bound on B {B ^ Kfi~^fi~), are not compatible with 
the experimental constraint on the B K*^^fi~ branching fraction. This is illustrated in 
Fig. ^, where we show the allowed range of -R9,io. Table |I| contains our predictions for the 
maximum average FB asymmetry and the branching ratios of -B — >■ K^*^ii'^n~ for certain 
choices of parameters. Note that the measurement of a nominal asymmetry of 4% (at 3a 
level), which is accompanied by a branching fraction of ~ 6 x IQ-'', will necessitate at least 
~ 10^'' B mesons and could conceivably be measured in forthcoming experiments at the 
CERN Large Hadron Collider (LHC) and the Tevatron. We conclude that the predicted FB 
asymmetry due to scalar interactions, though larger than in the SM, may possibly be too 
small to be seen experimentally. Nevertheless, the FB asymmetry does provide a very useful 
laboratory for studying possible extensions of the SM. 
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FIG. 4. The average FB asymmetry in S — > Kii^ijl~ as a function of (i?9,i?io) for Rs = —4, 
i?P = 0, and i?7 = 1 consistent with the upper hmit on B{Bs jjL^jjr). Also shown is the 
corresponding branching ratio B ^ Kn'^n" (see the text for details). 




-2-1012 -2-1012 

i?io Rio 



FIG. 5. Allowed ranges of Rg and Rio as determined from the upper limit onB{B ^ K*ii^n^) 
for Rs = —4, Rp = 0, and R7 = ±1. Note that the values of Rs,p are consistent with experimental 
data on Bg — ^ IJ-^ ■ 
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TABLE II. 


Maximum values of the 


average FB asymmetry in B - 




fi decay together 


with the branching 


ratios of B ^ 




for different choices of {R-j 


Rq,Rio). We have chosen 


{Rs,Rp) = {- 


4, 0), resulting in B [Bg — 




^fi~) ^ 2.5 X 10~^, which is 


close to the present upper 


bound of 2.6 x 10 


^. Also listed are the 


90' 


Vo C.L. upper limits as discussed in 


Sec. f^. 

I_J 


{Rt,Rq,Rio) 








B{B ^ Kfi+fi-) 




B{B ^ K*fi+fi-) 


(1,1.9,1) 




-2.6% 




1.5 X 10-6 




3.8 X 10-6 


(-1,1.2,1.3) 




-2.3% 




1.3 X IQ-^ 




3.9 X 10-6 


(1,1,1) 




-2.5% 




8.3 X 10"^ 




1.7 X 10-6 


(1,-1,0) 




+3.8% 




5.8 X 10"^ 




1.4 X 10-6 


(1,1,0) 




-3.8% 




5.5 X 10-^ 




7.8 X 10-^ 


(-1.2,1.1,0.8) 




-3.0% 




9.7 X 10^^ 




2.8 X 10-6 


Exptl limits 








< 5.2 X IQ-^ 




< 4.0 X 10-6 



B. Constraints on new physics with minimal flavour violation 

It is clear from the previous analysis that the upper bound on the Bs fi^^^ branching 
fraction gives the strongest constraints on the scalar and pseudoscalar couplings, cs,p, which 
in turn can be used to restrict the parameter space of models outside the SM. 



1. Two-Higgs-doublet model 

Using the constraints of the preceding sections, together with the results for the Wilson 
coefficients of scalar and pseudoscalar operators [Eq. (|5.9|) ], we obtain exclusions in the 
(m//±, Rs) plane, shown in Fig. ^. The observation to be noted here concerns the measured 
b ^ s'-f branching fraction, which gives the strongest constraint on the Higgs-boson mass, 
thereby placing an upper bound on \Rs\ of 0.17. As a consequence, the FB asymmetry due 
to non-standard scalar interactions turns out to be exceedingly small, typically of the order 
of 0(10-^). It should also be remarked that the lower bound on the charged Higgs-boson 
mass, mH± ~ 260 GeV, obtained from the measured inclusive b s'-f fraction in the context 
of the type-II 2HDM, does not directly apply to supersymmetric extensions of the SM since 
the chargino amplitude may interfere destructively with the charged Higgs- and ly-boson 
contributions. Given a charged Higgs-boson mass oi mH± = 260 GeV, the branching ratio 
of Bs —>■ fJ'^fJ'" amounts to (1.4-4.8) x 10-^ for 40 ^ tan/? ^ 60, which should be compared 
with the SM result of (3.1±1.4) x 10-^ [Eq. ( [4.17 )]. The average FB asymmetry is estimated 



to lie in the interval —0.7 x 10-'^ ^ (^fb) ^ —1.6 x lO-'*, which is much too small to be 
detected. For the decays B Kfi^fi" and B K*fi~^n~, we predict branching ratios 
of 5.2 X 10-^ and 1.7 x 10-6 respectively (to be compared with the SM expectations of 
about 6 X 10-^ and 2 x 10-6). Note that these decays are largely unaffected by the charged 
Higgs-boson contributions to Rg,Rio, which are proportional to cot^/3, and hence small in 
the large tan /5 regime. 

Our conclusion is therefore that current experimental data on various rare B decays - 
apart from 6 —> 57 - do not provide any constraints on the parameter space in two-Higgs- 
doublet models of class II. Moreover, the predictions for the branching ratios of the B decay 
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FIG. 6. Rs vs the charged Higgs-boson mass in the 2HDM for large values of tan /3, taking into 
account the constraint from Bs fJ-^fJ-^ decay. We exhibit the cases of tan /3 = 40 and tan f3 = 60, 
and also show, for comparison, lower limits on the Higgs-boson mass obtained by direct searches at 
the CERN e~^e~ collider LEP |^3| (vertical line at mjj± = 78.5 GeV) and by the measured 6 — > 57 
rate, using the leading-order result for the branching fraction. 



modes under study are comparable to those of the SM. We next turn to the SUSY scenario. 



2. SUSY with minimal flavour violation 

As mentioned earlier, we do not consider any CP-violating effects, and consequently the 
SUSY parameters and mixing matrices discussed in the previous section can be taken to be 
real. We further assume the sneutrinos to be degenerate in mass so that R^, which enters 
the expression in Eq. (|5.11|) , reduces to the unit matrix. 



For the sake of illustration, we perform the numerical analysis for a light stop ti, with 
large mixing 6^, and charginos with large Higgsino components. We impose the lower bounds 
on the SUSY particle masses as compiled by the Particle Data Group In the case of 



a light scalar top quark, there are additional constraints coming from electroweak measure- 



ments such as the p parameter ||3^. As for the constraint from 6 —> 57, it is well known that 
within supersymmetry there are chargino-stop contributions, in addition to charged Higgs 
boson and W boson loop contributions, which can significantly affect the b s'j decay 
rate in the large tan f3 region, thereby leading to constraints on the parameter space.0 This 



For example, within the minimal supergravity model, most of the parameter space is ruled out 
for /i < 0, where denotes the Higgsino mass parameter, while for /x > 0, much of the parameter 
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FIG. 7. Allowed region in the {mjj±,m-±) plane derived from the upper bound on 
_ ^1 
B (Bs fi^fi^) for a light stop with = 120 GeV and maximal mixing, 6^ ~ —45°. 



is due to the fact that at large tan/5, the chargino loop contributions grow linearly with 
tan/? (see, e.g., Ref. [^). Thus, in order to satisfy the b sj bounds, we must require 
that chargino and charged Higgs- and W^-boson contributions interfere destructively, so that 
significant cancellations can occur. Note that in this case the sign of the Wilson coefficient 
df is opposite to the SM one. 

Using the expressions for the scalar and pseudoscalar Wilson coefficients, Eqs. ( ^.9|) , 
( p.ll| )-( [5.1tj| ), we obtain the allowed (mj|/±,m-±) region displayed in Fig. |^ for m^-^ = 

120 GeV and 6^ ^ —45°. As can be seen, the present upper limit on B {Bg —>■ fi^fj^~) already 
excludes a significant portion of the SUSY parameter space with charged Higgs-boson and 
chargino masses. Remembering that the SM prediction for the branching ratio is of order 
O(10~^), it is clear that within the context of SUSY, the Bg fi~^fi~ branching fraction can 
be increased by several orders of magnitude, due to the tan^ /3 enhancement of the countert- 
erm diagrams [Eq. (|5.13|) ]. Given a chargino mass of m-± = 297 GeV, the SUSY prediction 
for Rs,p, consistent with the upper bound on B {Bs — > fi^fi"), is shown in Fig. ^( a^ , as a 
function of the charged Higgs-boson mass. From this we infer that Rs,p are constrained to 
lie in the range —3 < Rs,p ^ 3, leading to an average FB asymmetry of less than 2%. Figure 
P(b) displays the dependency of the Bg —>■ branching ratio on the charged Higgs-boson 

mass for m^± = 297 GeV. Thus, the measurement of the Bg ^^fi~ branching fraction, 
together with the 6 — 57 bounds, provides a useful tool for constraining supersymmetric 
extensions of the SM. Finally, for the above parameter space point and mH± = 170 GeV, we 



space is allowed by experimental data on 6 — > 57 [27,35| 
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FIG. 8. SUSY predictions for (a) Rs,p and (b) the branching ratio of Bg z^^^A*" at large 
tan/3, as a function of the charged Higgs-boson mass, with nif^ = 120 GeV, Of. « —45°, and 

m-± = 297 GeV. We note that this choice of parameters satisfies all experimental bounds. 

•^1 



obtain B {B ^ Kfj,+fx~) = 1.0 x 10"^ and B {B ^ K*fj,+fx-) = 2.6 x 10"^ to be compared 
with the present upper limits of 5.2 x 10~^ and 4.0 x 10~^. 



VII. SUMMARY AND CONCLUSIONS 

We have carried out a study of exclusive B decays governed by the b — > sl~^l~ transition 
in extensions of the SM with minimal flavour violation and new scalar and pseudoscalar 
interactions, focusing on the dimuon final state, and taking account of existing experimental 
data on 6 — > S7 as well as the upper limits on Bg — > fJ^'^fJ^' and B K^*^ n'^fjT decays. We 
have restricted the discussion to the interesting case of large tan /5, which may compensate 
for the inevitable suppression of scalar and pseudoscalar couplings by the lepton mass of e 
or /X. Our main findings can be summarized as follows. 

We have presented in a model-independent manner expressions for the Bg l^l" branch- 
ing fraction and the differential decay spectrum of B — > Kl~^l~, together with the corre- 
sponding FB asymmetry, which is extremely tiny within the SM. In particular, we find 
that scalar and pseudoscalar interactions can, in principle, lead to striking effects in the 
decay distribution of i? — > Kn'^n~, while the branching ratio of i? — > K*fj,^fj,^ is essen- 
tially unaffected by the Higgs-boson contributions. We have demonstrated that once the 
constraint from Bg iJ.^fjr is taken into account, the effects of scalar and pseudoscalar 
couplings on the decay B Kfi^fi~ are much smaller. In view of the inherent uncertainty 
of the predictions for exclusive B decays due to the form factors, it seems extremely unlikely 
that a measurement of the decay spectrum alone can provide a clue to new physics with 
scalar and pseudoscalar interactions. We have also investigated the FB asymmetry of /x~ 
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in B ^ K^'^fi~ decay, which turns out to be at the level of a few per cent. Our analy- 
sis suggests that the observation of a nominal FB asymmetry of, say, 4%, together with a 
branching ratio of about 6 x 10"'', will be challenging but might be feasible at the LHC and 
the Tevatron. As more precise data on the Bg A^"*"/^" branching ratio are available, more 
stringent upper limits will be placed on the FB asymmetry due to scalar interactions. The 
essential conclusion of our model-independent analysis is that current experimental data on 
Bg — > /i"*"/^" decay already exclude large values of the Wilson coefficients cs and cp of scalar 
and pseudoscalar operators, so that striking effects are not likely to show up in the decay 
spectrum of -B — Kfi'^n~ and the corresponding FB asymmetry of the muon. Even so, the 
FB asymmetry provides an independent window to physics beyond the SM, especially to 
models with an extended Higgs sector, and its observation would be an unambiguous signal 
of new physics. 

In extensions of the SM with minimal flavour violation, we have calculated the Higgs- 
boson contributions to the Wilson coefficients of scalar and pseudoscalar operators, and 
investigated how the new-physics parameters are constrained by existing experimental data 
on rare B decays. Within the type-II 2HDM framework, where the Higgs sector corresponds 
to the one of the MSSM, we found no appreciable FB asymmetry or any large deviation from 
the SM prediction for the B K^*^fj,~^n~ branching fractions. As for the decay Bg — >■ 
the branching ratio turns out to be in the range (1.4-4.8) x 10~^ for 40 ^ tan/5 ^ 60 and 
mu± = 260 GeV, which is within the errors of the SM prediction of (3.1 ± 1.4) x 10~^. 
Ultimately, the smallness of the new-physics contributions is caused by the mass of the 
charged Higgs boson, which is strongly constrained by the measured 6 — 57 branching 
fraction (m//± > 260 GeV). We therefore conclude that within the type-II 2HDM there are 
no sizable new-physics effects on the B decay modes described above, apart from h sj. 
By contrast, within SUSY, the effects of chargino and neutral Higgs-boson contributions on 
the Bg —>■ fi'^fJ'~ branching fraction can be enormous while satisfying the b ^ s'j bounds. 
We have considered a SUSY scenario with a scalar top quark with large mixing and a mass 
much lighter than the scalar partners of the light quarks. We find that for a given set of 
parameters obeying the b —>■ sj constraint, the upper limit on B {Bg — > /i^/i~) severely 
constrains the masses of charginos and charged Higgs bosons. As a typical result, the lower 
bound mjj± > 170 GeV has been derived for m-± = 297 GeV, m^^ = 120 GeV, and 
6f: ~ —45°. The remaining observables are estimated to be (Afb) = 1.8% and B{B ^ 
fi~^fi ) = 1.0 (2.6) X 10 (for mH± = 170 GeV), the latter being close to the present 
upper limits. Clearly, the analysis of the decay Bg — > is complementary to the study 

of i? — > K^*^ix~^fi~ and b ^ sj decays, which leads to constraints on the remaining short- 
distance coefficients, c'f, c^, and Ciq. A combined analysis of these decay modes, therefore, 
provides a powerful tool to constrain physics transcending the SM. 

Note added in proof. As this paper was readied for publication, we received a paper by 
Huang et al. [0 that corrects the result presented for the 2HDM in Ref . . The result of 
Ref. coincides with that of the present paper [Eq. ( |5.19| )]. 



ACKNOWLEDGMENTS 

We would like to thank Gerhard Buchalla, Andrzej J. Buras, and Manuel Drees for useful 
discussions. One of us (F.K.) would also like to thank Gudrun Hiller for communications. 



24 



We are indebted to Andrzej J. Buras for his comments on the manuscript. This work was 
supported in part by the German 'Bundesministerium fiir Bildung und Forschung' under 
contract 05HT9WOA0 and by the 'Deutsche Forschungsgemeinschaft' (DFG) under grant 
Bu.706/1-1. 

APPENDIX A: USEFUL FUNCTIONS 



1. The one-loop function Y{s) 

The function Y{s) appearing in Eq. ([4.13|) is given by |T4 



Y{s) = g{mc, s)(3ci + C2 + Scs + C4 + Scs + ce) - -g{0, s)(c3 + 804) 



where (at fib = m 



1 2 

-g{mb, s)(4c3 + 4c4 + Scg + Ce) + -(Scg + C4 + Scg + Ce) 

pole\ 



g 8 4 2 



X e(i-y,) 



In 



1 + 



Q{yi — 1)2 arctan 



with i/j = Amf/s, and 

ci = -0.249, C2 = 1.107, C3 = 0.011, C4 = -0.025, C5 = 0.007, Cq = -0.031. 



(Al) 



(A2) 



(A3) 



In Eq. (|A1|) , we have omitted the one-gluon correction to the matrix element of the operator 
Og, which can be regarded as a contribution to the form factors. 



2. Auxiliary functions 



B{x,y) 



C{x,y) 



y 



x-y 



Iny In X 



y 



Di{x,y,z) 



x-y 



xlnx 



y — I X — 1 



X In X y In y 
X — 1 y — 1 



^1 — x){x — y){x — z) 



D2{x,y) 



xlnx 



(1 - x)(x - y) 



+ [x ^ y) + {x ^ z) 



+ (x ^ y), 



D.ix) 



X Inx 
1 — X 



(A4) 
(A5) 
(A6) 
(A7) 
(A8) 
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APPENDIX B: SUSY MASS AND MIXING MATRICES 



For the reader's convenience and to fix our notation, we give here the relevant mass and 



mixing matrices in the context of SUSY, for which we adopt the conventions of Ref. p9[| 



1. Chargino mass matrix 

Neglecting new CP-violating phases, the chargino mass matrix is given by 
M^-( ^2 \/2MH.sin/3\ 

where M2 and fi are the VT-ino and Higgsino mass parameters respectively. This matrix can 
be cast in diagonal form by means of a biorthogonal transformation 

UM^±V^ = diag(m-±,m-±), (B2) 
m-± being the chargino masses with < m^±. The orthogonal matrices U, V read 

^1,2 00 Xi X2 

U = ( ^^^^^ V V=( '""'^l ^^^^7 1 , (B3) 

\ — sm Uu cos Uu I \ — sm Uy cos Uy I 



with the mixing angles 



2Mw[M^ + /i^ ± (M| - /i^) cos 2(3 + 2/iM2 sin 



1/2 



sin 2euy = " ' " " % — ' ^ ^g^^ 



— 2 2" 

171+ — m • 

Xi X2 



0. M|-/x^t2M^cos2/3 

C0S26^;7V' = ^ T • 

mt+ — mi+ 

Xi X2 



2. Scalar top mass matrix 

In the (t/i, J7ii) basis, the 6x6 up-type squark mass-squared matrix is given by 

f? M- 

Ulr Urr 

which can be diagonalized by an orthogonal matrix Rq such that 



Mi = J" J- I . (B6) 



flfjM^fli = diag(m4,ml,, . . . , mj,). (B7) 



Since we ignore flavour-mixing effects among squarks, the matrix in Eq. ( P6| ) decomposes 
into a series of 2 x 2 matrices. Working in the so-called super-CKM basis , in which the 



quark mass matrices are diagonal, and squarks as well as quarks are rotated simultaneously. 
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the LR terms in Eq. ( p6|) are proportional to ruq, q = u,c, t. Thus, large mixing can occur 
only in the scalar top quark sector, leading to a mass eigenstate m^^ possibly much lighter 
than the remaining squarks. 
Defining the 6x3 matrices 

(r^^),, = {R^),„ (r^«),, = {Ro)a,+3, (B8) 

we obtain 

10000 0\ /O 00100 

(r^^)^ = I 1 , (r^«)^ =00 01 | , (B9) 

cos^i -sinOi J \0 sin^^ cos^j 

with the mixing angle {—7t/2 ^ 6^ ^ 7r/2) 

mi — mi mi — mi 

Here At is the trilinear coupling, m^^ ^ are the soft SUSY-breaking scalar masses, and m^-^ ^ 
denote the stop masses with m? < m? . 

tl t2 

3. Sneutrino mixing matrix 

The 3x3 mixing matrix Rj^ appearing in Eq. (|5.11| ) is defined via 

R,M^Rj = dia.g{ml,ml,ml), (Bll) 



where M| is the sneutrino mass-squared matrix (see, e.g., Refs. ||38|,|39|). 
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